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Abstract 



The space-time evolution of open heavy flavor is studied in Pb+Pb collisions at y/s=2.76 TeV using the partonic transport 
jmodel Boltzmann Approach to MultiParton Scatterings (BAMPS). An updated version of BAMPS is presented which 
allows interactions among all partons: gluons, light quarks and heavy quarks. Heavy quarks, in particular, interact 
jwith the rest of the medium via binary scatterings with a running coupling and a Debye screening which is matched 
, * by comparing to hard thermal loop calculations. The lack of radiative processes in the heavy flavor sector is accounted 
O Tor by scaling the binary cross section with a phenomenological factor K = 3.5, which describes well the elliptic flow v-i 
and nuclear modification factor Raa at RHIC. Within this framework we calculate in a comprehensive study the V2 and 
'Raa of all interesting open heavy flavor particles at the LHC: electrons, muons, D mesons, and non-prompt J/ip from 
B mesons. We compare to experimental data, where it is already available, or make predictions. To do this accurately 
next-to-leading order initial heavy quark distributions are employed which agree well with proton-proton data of heavy 
' (—1 ' flavor at ^/s=7 TeV. 
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1. Introduction 

In ultra-relativistic heavy ion collisions at the BNL 
RHIC and CERN LHC the energy deposited in the col- 
lision zone is large enough to produce a medium that con- 
sists of deconfined quarks and gluons This state of 
matter, the quark gluon plasma (QGP), has remarkable 
properties such as collective behavior, a small viscosity to 
entropy ratio, and a large density that leads to quenching 
of jets. 

Open heavy flavor particles such as D and B mesons, 
which consists of one heavy quark and one light quark, 
are an exciting probe to study the QGP. Their heavy con- 
stituents, namely charm and bottom quarks, are created 
at a very early stage of the heavy ion collision due to their 
large mass [J]. Consequently, they travel for a long time 
through the QGP, collide with other medium particles, lose 
energy, and participate in the collective behavior. 

Measurements of heavy flavor electrons from the decay 
of D and B mesons at RHIC indicate that the elliptic 
flow V2 and nuclear modification factor Raa of open heavy 
flavor is on the same order as for light particles. This 
is in contrast to the expectations drawn from the dead 
cone effect [1, Q , that gluon radiation off heavy quarks is 
suppressed at small angles compared to light quarks. The 
reason for the rather large elliptic flow and suppression of 
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heavy flavor is currently under investigation 

At the LHC, for the first time, it is possible to dis- 
tinguish between charm and bottom quarks. In addition 
to looking at heavy flavor electrons or muons from both 



D and B mesons, ALICE can reconstruct D mesons di- 
rectly [23|]. In addition, CMS presented data for non- 
prompt J/ip [24[ which stem from the decay of B mesons. 

In the present paper we present our calculations of 
the elliptic flow and nuclear modification factor of D 
mesons and non-prompt J/ip, obtained with the trans- 
port model Boltzmann Approach to MultiParton Scatter- 
ings (BAMPS). In addition, we show results of heavy fla- 
vor muons at forward rapidity as well as electrons at mid- 
rapidity and compare them to experimental data. 

2. Open heavy flavor in BAMPS 

In this section we briefly review the features of our 
model, the parton cascade Boltzmann Approach to Mul- 
tiParton Scatterings (BAMPS). More details concerning 
the model itself and the implementation of heavy flavor 
can be found in Ref . (2|| [26| and [1, [2l[ , respectively. 

BAMPS is a 3+1 dimensional partonic transport model 
that solves the Boltzmann equation 
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for on-shell partons. Implemented processes on the light 
parton sector are all 2 — !► 2 and 2 o 3 processes. In con- 
trast to previous publications Jil], [27l429j | where we only 
took gluons (g) and heavy quarks (Q) into account, in the 
present calculation light quarks (q) are explicitly included. 
All cross sections are calculated in leading order pQCD. 



Preprint submitted to Elsevier 



March 4, 2013 



Light partons interact among each other via binary colli- 
sions and radiative 2 3 processes, which are taken in 
the Gunion-Bertsch limit [3Cj . 

For heavy quarks, currently only the elastic collisions 



gg^QQ 
qq o QQ 
gQ -> gQ 
qQ -> qQ 



(2) 



are implemented. The inclusion of radiative processes is 
underway and planned for the near future. 

In this paper we focus on the heavy flavor sector. For 
BAMPS results of light partons we refer to Refs. 2j| 26 
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The cross sections for the processes from @ are calcu- 
lated in leading order pQCD for a finite heavy quark mass 
37] . Since the matrix elements of the t channel of those 
elastic heavy quark scattering with a light parton are di- 
vergent, they are screened with a screening mass, which is 
determined from comparison to hard thermal loop (HTL) 
calculations. By comparing the energy loss of a heavy 
quark in a static medium calculated within HTL and the 
same quantity calculated from the leading order pQCD 
cross section with a screening mass fj, 2 = Km 2 n , one can 
obtain the prefactor k analytically to be 19|, 21 
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The Debye mass m 2 D is calculated in BAMPS from the 
non-equilibrium distribution functions / of gluons and 
light quarks via 25 1 
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where N c = 3 denotes the number of colors, v g = 16 is the 
gluon degeneracy, and 71/ = 3 the number of active light 
flavors. As a note, in equilibrium and with Boltzmann 
statistics the Debye mass is given by mp,^ = ^f-{N c + 
nf)T 2 . The Debye mass prefactor k was determined for 
an arbitrary number of light quark degrees of freedom n / 
and is thus easily applied to the t channel of heavy quark 
interactions with light quarks. 

Furthermore, instead of just assuming a constant value 
for the coupling a s we employ the running coupling for all 
heavy flavor processes [ijj 21. 38l 39$. 
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with fa = 11 - §71/ and L± = ln{±Q 2 /A 2 ) with A = 
200 MeV. For consistency a running coupling is also used 
in calculating the Debye mass from Eq. 

After the energy density in the surrounding of a heavy 
quark in BAMPS has dropped below 0.6GeV/fm 3 it is 



fragmented to a D or B meson. This is done according to 
the Peterson fragmentation function [4p| ] 
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N is a normalization constant, z = \ph\/\pq\ the ratio 
of the meson and quark momenta, and eg = 0.05 (0.005) 
for charm (bottom) quarks. The D mesons can then di- 
rectly be compared to the experimental data. To yield 
non-prompt Jh j) we carry out the decay of B mesons with 
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I42j by switching on the relevant decay chan- 
nels, pythia is also used to perform the decay of D and 
B mesons to electrons and muons which can then be com- 
pared to experimental data. 

Especially for the LHC, where charm and bottom can 
be separated, it is important to have the correct refer- 
ence for the initial heavy quark distribution. To generate 
the initial heavy quark spectrum for the BAMPS simu- 
lation of the heavy ion collision we employ the next-to- 
leading order (NLO) event generator MC@NLO [43], \4A \. 
The factorization and renormalization scales, fi F and /j,r, 
respectively, are in principle arbitrary when considering 
all orders of the cross section. However, for the lead- 
ing order cross section uncertainties due to neglecting 
higher order terms can be reduced if the two scales are 
of the order of the relevant scale \/p\ + M 2 , pp being 
the transverse momentum and M the mass of the pro- 
duced heavy quarks. The exact value of the scale is fixed 
by giving a good agreement with the experimental data 
which results in /ip = — 1 \/pp + M 2 for charm 
(Af c = 1.3 GcV) and [i F = \x R = 0.4 y/pj + M 2 for bot- 
tom quarks (M b = 4.6 GeV). 

In Fig. [T] the invariant differential cross sections of D 
mesons, heavy flavor electrons and muons are compared 
to experimental data from ALICE at y/s = 7TeV. The 
D mesons and heavy flavor electrons at mid-rapidity are 
well described by mc@nlo. At forward rapidity, however, 
the slope of the muons at larger pp is slightly different. 
Such a disagreement has also been observed by CMS in a 



48 



more detailed study of inclusive bottom jets in Ref. 
by comparing mc@nlo to data for larger pp and various 
rapidities. Nevertheless, we checked that the muon Raa 
is not very sensitive to the exact slope in this pp range. 

To obtain the initial heavy quark distribution as an in- 
put for BAMPS we run MC@NLO with the same parame- 
ters for a center of mass energy of </s = 2.76 TeV which 
was used for the recent heavy ion runs. 



3. Results 

Important observables for open heavy flavor are the nu- 
clear modification factor Raa and the elliptic flow V2- The 
nuclear modification factor is defined as the heavy flavor 
yield in heavy ion collisions divided by the yield from p+p 
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Figure 1: Differential invariant cross section of D mesons with \y\ < 
0.5 and heavy flavor electrons with \y\ < 0.8 at mid-rapidity and 
muons at forward rapidity 2.5<?/<4asa function of transverse 
momentum for proton-proton collisions with y/s = 7TeV simulated 
with MC@NLO. For comparison experimental data 45-47] with the 
same kinematic cuts is also shown. In the upper plot the electron 
curves and the corresponding data points have been scaled with the 
factor 0.1 to distinguish them from the muon curves. Since the data 
of electrons is preliminary, we do not have access to the errors and 
plot those data points without any errors as obtained from Ref . [4p| . 



collisions scaled with the number of binary collisions, 
d 2 N AA /dp T dy 



Raa 



N Un d 2 N pp /dp T dy 
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The elliptic flow denotes the second harmonic of the 
Fourier decomposition of the azimuthal particle spectrum 
and is given by 
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pI-pI 



For this definition the momenta p, 
respect to the reaction plane 
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and p y are taken with 
:8l we showed that elas- 



In previous publications [21I.I27I. 
tic heavy quark scatterings alone cannot describe the ex- 
perimental data at RHIC, although they play a significant 
role. To be compatible with the heavy flavor electron data 
for Raa and V2 at RHIC the elastic cross section had to be 
multiplied with a factor K = 4. Those calculations have 
been done without light quarks, nf = + 2 (number of 
light quarks + heavy quarks) . Figures [5] and [3] compare 
these calculations to the updated version of BAMPS which 
also includes light quarks {nf = 3 + 2). Comparing the 
K = 4 curves for nf =0 + 2 and u/ = 3 + 2 shows that 
both the suppression and elliptic flow are slightly higher for 
the latter. At first sight this seems counterintuitive since 
we use the same initial conditions from pythia for both 
cases and, to get the same energy density, convert initial 
light quarks to gluons for nf — + 2, which are associated 
with a larger Casimir factor. However, the situation is 
more complex. The running coupling also depends on the 
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Figure 2: Elliptic flow vi of heavy flavor electrons at RHIC with 
an impact parameter b = 8.2 fm together with data [7(. For heavy 
quarks only binary collisions are switched on, which are multiplied 
with a factor K. 
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Figure 3: Nuclear modification factor Raa of heavy flavor electrons 
at RHIC for the same configuration as in Fig. [2] 



number of flavors and is larger for n/ = 3 + 2. Further- 
more, the chemistry of a purely gluonic plasma behaves 
slightly differently as a quark gluon plasma. This leads to 
a larger number of scattering centers of the medium for 
nf = 3 + 2. All influences together result in a slight in- 
crease of the suppression and elliptic flow for nf = 3 + 2 
compared to nf = + 2. 

As can be seen in Figs. [5] and [3] the best agreement with 
data for nf = 3 + 2 is found with K — 3.5 in contrast to 
nf = + 2, where K = 4 yielded the best results. We as- 
sume that this factor is necessary due to the lack of radia- 
tive processes and quantum statistics in our calculations. 
It is planned to extend BAMPS in a forthcoming study 
to include also radiative processes. This will show if those 
processes can indeed account for such a phenomenological 
scaling of the binary cross section. As a notej+he value of 
3.5 is close to the needed K factor of Ref. [49(, which, in 
an independent framework, included similar cross sections 
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Figure 4: Nuclear modification factor Raa of heavy flavor electrons 
at Pb+Pb collisions at LHC with an impact parameter b = 3.6 fm 
together with data [5ll . For heavy quarks only binary collisions are 
switched on, which are multiplied with K = 3.5. 

for rif = 3 + 2 using ideal hydro as well as temperature and 
flow information from BAMPS for the medium evolution. 

For small px the employed hadronization scheme, 
namely Peterson fragmentation, is not valid and coales- 
cence might be the dominant process. It is expected that 
coalescence increases the elliptic flow at small transverse 
momenta, since light quarks would also contribute to the 
flow of the heavy mesons. This could be an explanation 
why BAMPS underestimates the flow in Fig. [5] for very 
small px- Neglecting cold nuclear matter effects such as 
the Cronin effect or shadowing and also coalescence effects 
could be the reason for the deviation of the Raa for small 
p T in Fig. [3] 

In conclusion, the effective description of the RHIC data 
with K = 3.5 agrees simultaneously with the data for 
both Raa and vi for intermediate and large px- In this 
paper we will employ this prescription at LHC energy of 
\fs — 2.76 TeV and compare to experimental data for D 
meson Raa and V2, non-prompt J/tp, electron and muon 
Raa- Furthermore, we make predictions for non-prompt 
J/ip, electron, and muon 1*2 . As a note, for all calcula- 
tions the same kinematic acceptance cuts are set at which 
experimental data is measured (see labels in the plots for 
the values). The impact parameters used in BAMPS are 
matched with a Glauber calculation to the mean number 
of participants (A part ) given for each centrality class (Hoj . 

In Ref. [2l[ we presented predictions for the electron 
Raa and at LHC for rif = + 2. In Fig. [4] an update 
of the heavy flavor electron Raa to nf = 3 + 2 is shown, 
employing again, as at RHIC, only collisional energy loss 
with a running coupling, improved Debye screening, and 
K = 3.5. Our calculation is consistent with the prelimi- 
nary experimental data, although it is on the lower edge 
of the uncertainty band. However, due to the rather large 
error bars one cannot judge yet whether the effective de- 
scription of the RHIC data also applies at the LHC. 



Figure 5: Raa of muons at forward rapidity at LHC with data 1531 . 

One huge advantage of transport models is the direct 
access to all particles of the QGP during the whole time 
evolution. Therefore, we plot in addition to the heavy 
flavor electron curve the Raa of charm and bottom quarks. 
As can be seen in Fig. 2] the curve of the heavy flavor 
electrons is shifted to smaller pt compared to the parental 
heavy quarks due to fragmentation and decay processes. 
Predictions for electron elliptic flow will be presented along 
with muons and non-prompt J /ip at the end of the section 
since there is no experimental data available yet. 

ALICE can measure muons in forward rapidity stem- 
ming from heavy flavor decays. Figure [5] shows the 
BAMPS results in the same rapidity range as the experi- 
mental data. Comparing to Fig. |4j it is obvious that the 
suppression of muons at forward rapidity is as strong as 
that of electrons at mid-rapidity and both nuclear mod- 
ification factors assume very similar values. In contrast 
to the electron data, the muon data has only small errors 
and a deviation is visible between the data and our curve, 
which is calculated for the same parameters that describe 
the RHIC heavy flavor electron data. In addition to the 
muon curve, we show the Raa on the heavy quark level 
directly which is again shifted to larger px as one would 
expect. 

By considering muons or electrons, the contributions 
from charm and bottom cannot be distinguished. How- 
ever, at the LHC for the first time one has access to charm 
and bottom separately via D mesons and non-prompt J/ip 
from B mesons. 

Figure [6] depicts the BAMPS calculations for D meson 
V2 ■ The error bars of the preliminary data are too large to 
draw any definite conclusion, but our results are in good 
agreement within the errors. In addition to the D meson 
curve we plot the charm quark 1*2 . The difference between 
both curves is considerably smaller than it was the case for 
heavy flavor electrons or muons, which renders D mesons 
a very good indicator for actual charm quark observables. 

The D meson Raa from BAMPS is compared to data in 
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Figure 6: Elliptic flow i>2 of D mesons at Pb+Pb collisions at LHC 
with an impact parameter b = 9.7 fm together with data [53|. 



Figure 8: Raa of non-prompt J/jf> at LHC with data [24 
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Figure 7: Nuclear modification factor Raa of D mesons at LHC with 
data [H. 



Fig. [7] Although the order of magnitude of the suppression 
is comparable, the experimental data tends to be slightly 
underestimated by our calculation. This is in accordance 
to the muon Raa at forward rapidity in Fig. [5j which is 
also below the data, and the electron Raa in Fig. 01 which 
is at the lower edge of the error bars. This could be a first 
hint that new effects compared to RHIC play a role at the 
LHC. An indication in this direction is also the fact that D 
meson suppression seems to be slightly smaller than that 
of charged hadrons (23j . 

Possible explanations for the discrepancy in our Raa 
calculations and the heavy flavor data could be cold nu- 
clear matter effects, the normalization error of the data 
which is not shown in the plot or that we represent the 
rather large centrality classes by only one impact param- 
eter. Furthermore, a reason could be that the approxima- 
tion of modeling the radiative energy loss by scaling the 
binary cross section with a constant factor is not satisfied. 
Although we do not expect that such a K factor is tem- 



perature dependent for a thermalized system, non-thermal 
effects in the medium evolution could trigger different K 
factors for different collision energies at RHIC and LHC. 
However, this cannot be assessed without actually doing 
the calculation with higher order processes where the K 
factor is obsolete. We will investigate this in more detail 
in a forthcoming study. 

A complimentary measurement has been performed by 
the CMS collaboration [2^ which measured the suppres- 
sion of non-prompt J/tp from the decay of B quarks. Al- 
though only one data point could be extracted, the sup- 
pression of non-prompt J/ip is clearly visible in Fig. [5] and 
the magnitude is in good agreement with our calculation. 
Analogously to the other Raa comparisons at LHC, our 
curve is slightly smaller than the experimental value, al- 
though still within the errors. Again, the suppression of 
bottom quarks themselves is very similar to that of non- 
prompt J/ip. 

To conclude the vi and Raa comparisons we show in 
Fig. IH1 BAMPS predictions of the elliptic flow of muons, 
electrons, and non-prompt J/ip calculated with the same 
parameters used for the previous figures, which describe 
the RHIC data. For better comparison the curve of D 
mesons from Fig. [5] is also depicted. 

The flow of non-prompt J/ip is considerably smaller 
than the D meson flow due to the mass difference of charm 
and bottom quarks. Accordingly, the influence of bottom 
quarks to the flow of heavy flavor electrons at intermediate 
and large pt is also the reason why the electron flow docs 
not increase to the value of the D meson flow. Muons 
at forward rapidity adopt the same elliptic flow as elec- 
trons at mid-rapidity. This is in accordance with the same 
Raa of muons and electrons (cf. Figs. [5]and@|. Since 
BAMPS is a 3+1 dimensional transport model, boost in- 
variance of the system in rapidity is not assumed, but 
in first approximation - comes out naturally for not too 
large rapidity gaps, which is reflected in the same vi and 
Raa of electrons and muons at mid- and forward rapidity, 
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Figure 9: Predictions for the elliptic flow vi of muons, electrons, 
non-prompt J/ip, and D mesons. 



respectively. 

Finally, the question arises, what we can learn from the 
comparison to experimental data. Due to the transport 
character of BAMPS we have access to all the collision 
properties during the whole time evolution. Figure QJJ 
sheds some light on why the experimental data can be 
fairly well described with the K factor. 

For the elliptic flow isotropization is important and, 
hence, the transport cross section and transport rate are 
the relevant quantities, since they weight the cross section 
and rate, respectively, with the angle of the diffracted par- 
ticle. In the left plot of Fig. [TUJ the time evolution of the 
mean transport cross section is shown in the central re- 
gion of a heavy ion collision at LHC. The value for charm 
quarks, including only 2 — > 2 processes, is about 5 times 
larger than that for gluons, which interact also via 2h3 
processes. As a note, the 2 — > 2 and 2 <H- 3 light parton 
processes can describe the elliptic flow data of light parti- 
cles 
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35| . Due to the large charm transport cross 
section also a sizeable elliptic flow for charm quarks builds 
up. 

For comparison we show also the charm transport cross 
section under the assumption that the medium is ther- 
mally and chemically equilibrated. These values are 
smaller compared to the full heavy ion collision, where 
the gluon fugacity is below unity, which leads to a smaller 
Debye mass and, therefore, larger cross section. However, 
the sensitivity on the fugacity is reduced when considering 
the transport rates since the density also enters. The right 
hand side of Fig. QJJ shows that for the charm transport 
rate the equilibrium curve is only a factor of about 1.4 
smaller than the values extracted from the full heavy ion 
collision. 

The reason for the large charm cross section is the ef- 
fective treatment of the Debye screening on the heavy fla- 
vor sector and the additional multiplication of the factor 
K = 3.5. Since the binary pQCD cross section is domi- 
nated by small angles, the total cross section can be one 



order larger than the transport cross section, which is too 
large for partonic interactions. Therefore, we will inves- 
tigate in a further study radiative processes, which can 
lead to an effective energy loss and isotropization (cf. also 
Ref. 22]) without the need of large cross sections, as it 
has been demonstrated with gluons and light quarks in 
BAMPS [H. 

4. Conclusions 

Results on the elliptic flow and the nuclear modifica- 
tion factor of open heavy flavor at the LHC are calcu- 
lated within the partonic transport model BAMPS with 
rif = 3 + 2 flavors. To get the correct initial heavy quark 
distribution we employ MC@nlo which is in good agree- 
ment with p+p heavy flavor data at y/s — 7 TeV. The 
binary cross sections of heavy quarks with light medium 
particles are obtained from pQCD with a running coupling 
and an improved Debye screening motivated from hard 
thermal loop calculations. By scaling the binary cross sec- 
tion by a factor of K = 3.5 to account for radiative contri- 
butions and quantum statistics we find a good agreement 
with the heavy flavor electron data at RHIC. Calculations 
with the same parameters are carried out at the LHC for 
heavy flavor electrons, muons, D mesons and non-prompt 
J/ip. We find that our Raa calculations slightly underes- 
timate the experimental data for all heavy flavor observ- 
ables. In a future study we want to investigate this in more 
detail by explicitly including radiative processes for heavy 
quarks [9( and study whether a phenomenological K factor 
was satisfied. The D meson t>2 agrees well with the exper- 
imental data. Furthermore, we make predictions within 
the same framework for the elliptic flow of heavy flavor 
electrons, muons, and non-prompt J/ip at LHC energy of 
s/s = 2.76 TeV. 
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all particles in a tube with a radius of 1.5 fm and space-time rapidity rj £ [—0.5; 0.5]. In addition, the charm transport cross section and rate 
in a chemically equilibrated medium is shown. That is, for each time we extract the medium energy density and mean charm energy from 
BAMPS in the central region and compute the transport cross section and rate of a charm quark with that energy in a static and equilibrated 
medium. 
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